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Abstract. In this paper, the Choquet integral and the interval neutrosophic set theory are combined to make multi-criteria decision 
for problems under neutrosophic fuzzy environment. Firstly, a ranking index is proposed according to its geometrical structure, and 
an approach for comparing two interval neutrosophic numbers is given. Then, a < L implied operation-invariant total order which 
satisfies order-preserving condition is proposed. Secondly, an interval neutrosophic number Choquet integral (INNCI) operator 
is established and a detailed discussion on its aggregation properties is presented. In addition, the procedure of multi-criteria 
decision making based on INNCI operator is given. Finally, a practical example for selecting the third party logistics providers is 
provided to illustrate the feasibility of the developed approach. 
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1. Introduction 

The concept of neutrosophic set (NS) is, intro- 
duced by Smarandache [22], which generalizes the 
classic set, fuzzy set (FS), interval valued fuzzy set 
(IVFS), intuitionistic fuzzy set (IFS), as well as inter- 
val valued intuitionistic fuzzy set (IVIFS). A NS is 
characterized independently by a truth-membership, an 
indeterminacy-membership and a falsity-membership. 
It is a powerful tool to deal with incomplete, indeter- 
minate and inconsistent information. Comparing with 
NS, IFSs and IVIFSs can only handle incomplete 
information but not the indeterminate information and 
inconsistent information which exist commonly in real 
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situations. For example, in a decision making pro- 
cess, a manager decides whether he should select the 
third party logistics provider A or not. When we ask 
about the opinion of an expert about a certain state- 
ment, he may say that the possibility that he select A 
is between 0.5 and 0.6, that he does not select A is 
between 0.1 and 0.2, and the degree that he is not sure 
is between 0.2 and 0.3. For a neutrosophic notation, it 
can be expressed as ^([0.5, 0.6], [0.1, 0.2], [0.2, 0.3]). 
The sum of the degree of truth, indeterminacy, and 
falsity may be greater or less than 1. For example, 
if x — (0.6, 0.2, 0.3), x is a NS but is not an IFS. If 
x= (0.5, 0.1, 0.2), x is a NS but is not an IFS. If 
x = (0.5, 0.2, 0.3), x is not only a NS, but also an IFS. 
Another example, assuming there are 10 voters during 
a voting process, in time t\, four vote ’yes’, three vote 
’no’ and three are undecided. For neutrosophic notation, 
it can be expressed as a'(0.4, 0.3, 0.3); in time t%, two 
vote ’yes’, three vote ’no’, two give up, and three are 
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undecided, then it can be expressed as *(0.2, 0.3, 0.3). 
That is beyond the scope of the IFS. So the notion of 
neutrosophic set is more general [22, 33]. 

In aNS, the degree of truth, indeterminacy, and falsity 
belong to ]0“, l + [, where ]0 _ , l + [ is the non-standard 
unit interval [21]. Obviously, it is difficult to apply 
in real applications. Therefore, Wang et al. [27] pro- 
posed the concept of a single valued neutrosophic set 
(SVNS), which is the subclass of a NS. Sometimes the 
degree of truth, falsity, and indeterminacy of a certain 
statement cannot be defined exactly in the real situa- 
tions but denoted by several possible interval values. 
In order to research this problem, Wang et al. [26] pro- 
posed the concept of interval neutrosophic set (INS) and 
gave the set-theoretic operators of INS. Recently, many 
researchers have shown great interest in multi-criteria 
decision making (MCDM) problems with neutrosophic 
information. Ye [30, 31] proposed correlation coeffi- 
cients between SVNSs and applied them to MCDM 
problems with single valued neutrosophic informa- 
tion. Ye [32] proposed single valued neutrosophic cross 
entropy and applied it to MCDM. Ye [29] introduced the 
concept of simplified neutrosophic sets (SNSs), which 
can be described by singleton subintervals/subsets in 
the real unit interval [0, 1], and proposed a MCDM 
method using aggregation operators for SNSs. Chi and 
Liu [2] extended a TOPSIS method to interval neu- 
trosophic multiple attribute decision-making problems. 
Zhang et al. [34] defined the operations of INSs and 
gave the aggregation operators of interval neutrosophic 
number weighted averaging(INNWA) and interval neu- 
trosophic number weighted geometric(INNWG), then 
a MADM method is established based on the proposed 
operators. Broumi and Smarandache [1] introduced the 
concept of correlation coefficients of interval valued 
neutrosophic set. 

Aggregation function plays an important role in 
MCDM problems. All above aggregation operators 
only consider situations where criteria (attribute) and 
preferences of decision makers are independent of one 
another, which means that their effects are viewed as 
additive. However, in many real decision making prob- 
lems, it is common to find that there is interaction among 
preference of decision makers. As an aggregation func- 
tion, the Choquet integral [3] with respect to fuzzy 
measures [23] is able to flexibly describe the relative 
importance of decision criteria as well as their interac- 
tions [11, 19]. Therefore, it is interesting to combine 
the Choquet integral and the INS theory for MCDM 
under neutrosophic fuzzy environment. On one hand, 
we can deal with the imprecise and uncertain decision 



information; on the other hand, we can efficiently take 
into account of the various interactions among the deci- 
sion criteria. Based on the above discussion, there are 
three aims in the paper. First, it will propose a new rank- 
ing method for interval neutrosophic numbers (INNs). 
Second, it will propose an interval neutrosophic num- 
ber Choquet integral (INNCI) operator and discuss its 
properties. Third, it will establish a decision making 
method based on the proposed ranking method and the 
INNCI operator to handle decision making problems 
with interval neutrosophic information. 

The paper is organized as follows. In Section 2, the 
concepts of NS, SNS, INS and operations of INS are 
reviewed. In Section 3, a ranking index is proposed 
according to geometrical quantities which reflect the 
inferiors and superiors of INNs, and an approach to 
compare two INNs is proposed. Furthermore, a <i 
implied operation-invariant total order which satisfies 
order-preserving condition is proposed. In Section 4, 
INNCI operator is proposed, and some of its proper- 
ties are researched. In Section 5, the MCDM procedure 
based on INNCI operator is presented under neutro- 
sophic environment. In Section 6, an example is given 
to illustrate the concrete application of the method and 
to demonstrate its feasibility and applicability. Conclu- 
sions are made in Section 7. 



2. Interval neutrosophic Set 

Definition 2.1. [22] Let X be a space of points ( objects), 
with a generic element in X denoted by x. A NS A 
in X is characterized by a truth-membership function 
Ta(x), an indeterminacy-membership function Ia(x), 
and a falsity -membership function Fa(x). The func- 
tions Ta(x), I a(x) and Fa(x) are real standard or 
nonstandard subsets of ]0 _ , l + [, i.e., Ta(x) : X — > 
]0“, l + [, Ia(x ) : X ->]0“, l + [, and Fa(x) : X — »■ 
]0“, l + [, where ]0 _ , l + [ is the non-standard unit 
interval. The sum ofTA(x), Ia(x ) and Fa(x ) satisfies 
the condition 0“ < sup Ta(x) + sup /a O r) + sup Fa(x) 
< 3+ 

Since it is difficult to apply NSs to practical problems, 
Wang et al. [27] introduced the concept of a SVNS, 
which is an instance of a NS and can be used in real 
scientific and engineering applications. 

Definition 2.2. Let X be a space of points ( objects) with 
generic elements in X denoted by x. A SVNS A in X 
is characterized by truth-membership function Ta(x), 
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indeterminacy-membership function I a (x), citid falsity- 
membership function F A (x). For each point x in X, 
Ta(x), Ia(x), Fa (x) e [0, 1]. A SVNS A can be written 
as A = { (x, T a (x), Ia(x), F a (x))\x e X}. 

Similar to interval-valued intuitionistic fuzzy set, 
Wang et al. [26] proposed the concept of INS. 

Definition 2.3. [26] Let X be a space of points ( objects) 
with generic elements in X denoted by x. An INS A in X 
is characterized by a truth-membership function Ta(x), 
an indeterminacy-membership function I A (x), and a 
falsity-membership function Fa(x). For each point x 
in X, we have that T A (x) = [inf T A (x), sup Ta (x)] C 
[0, 1], I a (x) = [inf / A (x), sup I A (x)] c [0. 1], F A (x) = 
[inf F a (x), sup F a (x)] c [0, 1], 

Remark 2.1. From Definition 2, an INS A can be 
expressed as 

A — {(x, 7a(x), I a (x), F a (x)) | x e X} 

= {{ x , [inf 7 a (x), sup T A (xj], [inf 7 a (x), sup ! a (x)]. 

[inf F a (x), sup F a (x)]} \ xe X} 

Then, the sum ofT A (x), I A (x) and F A (x) satisfies the 
condition 0 < sup 7 a (x) + sup 7 a (x) + sup F A (x ) < 3. 

Definition 2.4. [26] An INS A is contained in the 
other INS B, A C B, if and only if inf 7 a (x) < inf 
7g(x), sup 7a ( x) < sup Tb(x), inf 7a ( x) > inf 7g(x), 
sup 7a ( x) > sup I p(x), inf F a (x) > inf I-ffix), and 
sup F a (x) > sup Fg{x) for any x in X. 

For convenience, let a = ([7^', T~ J ], [iff, if], 
[Ff , fV]) denotes an INN and L be the set of all INNs 
in X. 



k(x) — — log(x), the operational rules defined in [2] and 
[34] are identifiable. 

Theorem 2.1. Let d and b be two INNs, then a 0 b and 
X <S> a are also INNs. 

The results of Theorem 2.1 are obvious. Further- 
more, the operations of two INNs have the following 
properties [34], 

Proposition 2.1. Let a ,b and c be INNs and a. f> > 0, 
then 

(1 )d@b = b@d,d®b = b®a. 

(2) (a © b) © c — a © (b © c), (d ® b) ® c = a ® 

( b (8) c). 

(3) a(d ®b) = ad(B ab, ( a © b) a = d a © b a . 

(4) ad © fa — (a + fi)a, d a <g> a^ = d a+ P . 

3. Comparison method of INNs 

In this section, we will propose two approaches to 
compare two INNs. Firstly, a ranking index is devel- 
oped according to some geometrical quantities which 
reflect the inferiors and superiors of INNs, and then 
order relation <# is proposed. Furthermore, due to lim- 
itations of <h, a <L implied operation-invariant order 
is proposed. 

3.1. Ranking method for INNs based on 
geometrical structure 

Let us consider a three-dimensional coordinate of 
the technical neutrosophic cube that shown in Fig. 1, 
where x is the truth axis with value range in [0, 1], y is 
the false axis with value range in [0, 1 ], and similarly z 



Definition 2.5. [2] Let d and b be two INNs, and X be 
a real number. Then, the operational rules are defined 
as follows: 

(1) d®b= ([t1- 



T b 



'j'L r rL r rU I r rU 

a ' 1 b ~ r b 

T U a ■ ], ut ■ I ! b X 1% • 7f], [Ft • F^, fV ■ fV\). 

(2) 5 ® b = ([T-t ■ T< b , TlJ ■ T-V], [It + I-t ~ I-t • 
jL jU , jU _ jU . tU i r pL , pL _ ri . pL pU , 
b ' a ' b a b *’ L a ' b a b ’ 5 ' 

pf _ pU . pUc, 
b a b u 



(3) Xd = ([1 - (1 - T-t)\ 1 - (1 - Tttl [dt)\ 
(7a / ) X n(^ i )\(7f) A ]),k>0. 

(4) = ([(r a L )\ (Tl/y i [1 - (1 - It)\ i — (i — 

I’in [i - (i - Ft)\ i - (i - F%f]), x > o. 



The operations of two INNs has been defined by 
Zhang et al. [34] by using of a strict Archimedean t- 
norm T(x, y ) = k~(k{x) + k(y)) and its dual t-conorm 
S(x, y) = l~(l(x) + l(y)) with Z(x) = k(l — x). When 




Fig. 1. The three-dimensional coordinates of the INNs. 
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Fig. 2. The two-dimensional projection of INNs. 



From the Definition 3 . 1 , it is easy to get the following 
proposition. 

Proposition 3.1. For any INNs a = ( \Tf, Tf ] , 
[if, if, [Ff, Ff] >, then H(a) e [0, 1]. Further- 
more, if a — ([1, 1], [0, 0], [0, 0]}, then H(d) — 1, and 
ifd=([ 0, 0], [1, 1], [1, 1]}, then H(d ) = 0. 

Theorem 3.1. Let a = ([T-f, if], [if. if], 
[Ft, Ft]) and b = ([7* T-f], [if if, [Ff Ff]) 
be two INNs. If T-t < T-f, Tf < T-f, if > if 
it > if, Ft > F-f and Ft > F- 7 , f/ten 77(a) < 77(b). 



is the indeterminate axis with value range in [0, 1]. The 
neutrosophic cube can be divided into three disjoint 
regions [5]: 

(1) Triangle BDE , whose sides are equal to V 2, 
represents the geometrical locus of the points whose 
sum of the coordinates is 1 . 

(2) The pyramid EADB is the locus of the points 
whose sum of coordinates is less than 1 . 

(3) the solid EHDCGFB ( excluding A BDE ) is the 
locus of points whose sum of their coordinates is greater 
than 1. 

Considering a point a in the technical neutrosophic 
cube. The superiors point is 5(1, 0, 0) and the inferi- 
ors point is HU). 1,1), therefore, the shorter distance 



Proof. Let G(d) = 



S i + D* 

s 2 + d- 




between a and 5(1, 0, 0) is, and the longer distance 
between a and H( 0, 1, 1) is, the bigger a is. 

Given a INN a, the two-dimensional projection of a 
in plan (x. y) is shown in Fig. 2. The smaller area of Si 



s 2 + d- 



Sj + D* 1 

(S 2 + D~f X 2 





is, and the bigger area of S 2 is, the bigger d is. The 
shorter distance between (rf + Tf /2, if + if /2, 
fF + Ff / 2) and 5(1, 0, 0) is, and the longer distance 
between ( T-f + T-f/2, if + if/ 2, Ff + F f /2 ) and 
7/(0, 1 , 1) is, the bigger a is. According to the geomet- 
rical description of INNs, the ranking index of INNs is 
proposed. 



dG(a) 

Since T a C [0, 1], F a C [0, 1], then - 

8G(a) 8G(a) 8G(a) 

0, — V < 0. Similarly, — -f > 0 , — \± > 0, 
8T- U ~ dlF ~ dlF ~ 

u a a u a 

^ n ^ n 

aA - ’ a vU - U ' 



Definition 3.1. Let d be an INN. The ranking index of 
a is defined as follow 



H(d ) = 



S 2 + D~ 

Si + S 2 + D* + D- 



where Si = (1 - T-f) ■ F-f , S 2 = (1 - 

Ff) ■ Tf , D* = yJ(l-Ttf + ll + Fl 

D- = \/t£ + (1 — la) 2 + (1 — Fa) 2 , Ta = 



'T'L I r rU 

F n ' I 7, 



— T~ — 

9 1 a — 



iF + iV 

n ' n 



— and Fa = 



FF + fF 

n 1 n 



Thus, G(d) is a decreasing function of Tf and if , 
an increasing function of if, if ) ff and Ff . There- 
fore, 1 1(d) is an increasing function of Tf and if , a 
decreasing function of if , if , Ff and Ff . The proof 
of Theorem is completed. 

Based on the ranking index of INNs, an approach to 
compare two INNs is proposed as follow. 

Definition 3.2. Let d and b be two INNs, Hid) and 
H(b) be ranking index of a and b, respectively, then 
(1) If H(d) < H(b), then d is smaller than b, denoted 
by a <h b- 
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(2) If H (a) — H(b), then a and b represent the same 
information, denoted by a —h b. 

The order relation <# is reflexive, antisymmetric, 
transitive and total, and hence defines a total order on 
INNs. 

3.2. </ implied operation-invariant order 



•IfT’/ < 



TP then a <g b; 
• If TP = T-f, then 



< T- 



♦ IfTt = f 

7 U 



IfF-a < F h 
If fV = F f 



TP 1, 

a J’ 



uk 



/'/]• [Ft 



Definition 3.3. Let a — ([Tf 
FP ]) andh= {[TP , TP], [if;, if], [Ff if]) be two 
INNs. An order relation <L on L is defined by 

( 1 ) a < L b iff TP < T-f, tP < T-f , / a L > if 



® If FP < Ft 
®IfFP = F-, 



then a <$ b; 
then 

u then d <s b; 
then 

then ci <s b; 
then 



*IfI~a <I- b > 

* * If iP < if, then a <s b; 



L 
b ’ 
L 
b ’ 
jU 



then a <s b; 
then 



if > if F P > F-P and Ffi > F~ 



r u 



U 



?u. 



b 



( 2 ) a < L b iff a < b and a fib 



The order relation <p is reflexive, antisymmetric and 
transitive, and so a partial order. As compared with 
the order relation <h, the order relation </ has some 
important properties although it is not a total order. 
According to Theorem 3.1 and Definition 3.3, it is easy 
to obtain the following proposition. 



b ’ 

★ If FP = FP then a =5 b. 

It is easy to verify that the order relation <s satis- 
fies the implication conditions and the order-preserving 
conditions given in the Definition 3.4. 

3.3. Comparative analysis with score, accuracy 
and certainty functions for INNs 



Proposition 3.2. Let a, b and c be INNs and a > 0. 

(1) If a <l b. then a <h b. 

(2) If ci <l b, then a@c <l b ® c and ci ® c <l 
b ® c. 

(3) If a <p b. then ad <l ab and d a <p b a . 
Monotonicity is one of the most important properties 

of an aggregation function. Therefore it is necessary to 
define an operation-invariant total order. Inspired by 
the concept of the operation-invariant total order on 
IFVs proposed by Liu [15-17] and an order implied 
operation-invariant total order proposed by Wu et al. 
[25], the following definition is given. 

Definition 3.4. Let a, b and c be INNs and a > 0. A 
<1 implied operation-invariant total order, denoted by 
<, is a total order on INNs if it satisfies the following 
conditions 

(1) If a <L b, then d<b. 

( 2 ) Ifa<b, then a © c<b © c and a ® c<b < 8 > c. 

(3) Ifa<b, then ad<ab and d a <b a . 

In the following, we will propose a <p implied 
operation-invariant total order. 



J 

"V 



Zhang et al. [34] defined the score function, accu- 
racy function and certainty function for an INN a — 
<[ TP , TP ], [ IP , If], [FP, fP ]) as follows: 



s(ci) = 



TP+ 1 



+ 1 -^ + 1 



i(d) = [min {TP - FP 
min{TP-FP, 
c(d) = [TP 



1 — fV tV 



T~ — F~ 



T~ — F~ 



J^U 



(2) 

(3) 

(4) 



Definition 3.5. Let a = ([if, TP], [if, if]. [FP, 
Ff ]) andb = ([Tf, T-f], [if, if], [F-f, Ff ]) be two 
INNs. A <1 implied operation-invariant total order, 
denoted by <g, can be defined as follows: 



where s(d), a(a) and c(a) represent the score function, 
accuracy function, and certainty function of the INN d, 
respectively. 

Zhang et al. [34] gave the ranking method as follows: 

Definition 3.6. Let a and b be two INNs, then 

(1) If p{s(d) > s(b)) > 0.5, then a is greater than b, 
denoted by d > p b; 

(2) If p(s(a) > s(b)) — 0.5 and p(a{d) > a(b )) > 
0.5, then a is greater than b, denoted by a > p b; 

(3) If p(s(a) > s(b)) = 0.5, p(a(d) > a(b)) = 0.5 
and p(c( a) > c(b)) > 0.5, then a is greater than b, 
denoted by d >p b; 

(4) If p{s{a) > s(b)) = 0.5, p(a(d) > a(b )) = 0.5 
and p(c(d) > c(b)) — 0.5, then d is equal to b, denoted 
by a = p b. 
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where P(a > b ) denotes the degree of possibility of a > 
b and is formulated for a = \a L , a u ] and b = [b L , b u ] 
by 



P(a > b) — max 



1 — max 



b u -a L 

a u — a L +b u — b L 





Remark 3.1. The order relations <h and <p can 
reflect the characteristics of INNs and have their own 
advantages. The order relation <h puts emphasis on 
geometrical structure of INNs, and <p can reflect 
the possibility degree of a > b. Comparing with <h 
and <p, the order relation <g puts a heavy emphasis 
on the truth-membership degree of an INN and more 
or less neglects the indeterminacy-membership degree 
and falsity-membership degree. It does not adequately 
reflect the characteristics of INNs. 

The order relation <h and <p have some limitations 
which will be shown in examples below. 

Example 3.1. Let a = ([0.65, 0.75], [0.10, 0.30], 
[0.25, 0.55]}, b = ([0.55, 0.85], [0.10, 0.30], [0.35, 
0.45]}. 

(1) The degrees of possibility of s(ci) and s(b )), a{a) 
and a(b)), c(a) and c(b )) are computed as follows: 

p(s(a) > s(b)) = p(a(a) > a(b)) — p{c(a) 

> c(h)) = 0.5. 

then from Definition 3.6, a =p b. But it is obvious that 
a and b are not equal. 

(2) By Eq. ( 1 ), H(a) = 0.7157, H(b ) = 0.7205, then 

a <H b. 

Example 3.2. Let a = ([0.65, 0.75], [0.20, 0.45], 
[0.32, 0.38]}, b = ([0.62, 0.68], [0.25, 0.40], [0.25, 
0.35]}. 

(1) The degrees of possibility of s(d) and s(b)), a{a) 
and a(b)), c(a) and c(b)) are computed as follows: 

p(s(d) > s(b)) — p(a(a) > a(b)) — 0.5, 

p(c(a) > c(b)) = 0.8125. 

then from Definition 3.6, a > p b. 

( 2) ByEq . (1), H(a ) = H(b) — 0.7095, then a —h b- 
But it is obvious that ci and b are not equal. 

The order relation <h and <p are not order- 
preserving for operations a © b, a <g> b, Xd and a x . In 
the following, counterexamples are given. 

Example 3.3. Let a = ([0.66, 0.71], [0.35, 0.26], 
[0.26, 0.30]}, b = ([0.66, 0.70], [0.32, 0.42], [0.20, 



0.23]}. From Definition 4.3, H(a) = 0.7351, H(b) = 
0.7435, then a <h b. From Definition 3.6, p(s(a ) > 
s(b)) — 0.4848, then ci <pb. 

(1 )Letc= ([0.64,0.90], [0.51,0.52], [0.10,0.70]}, 
then a@c>Hb@c,a@c>pb@c. 

(2 ) Letc = ([0.64, 0.90], [0.51, 0.52], [0.68, 0.70]}, 
then a®c>Hb<g)c, a<g>c>pb®c. 

(3) Let X — 2, then Xd >h Xb, Xd >p Xb. 

(4) Let X — 0.2, then a x >h b k , then a x >p b x . 

Remark 3.2. From Example 3.1 and Example 3.2, in 
some situations, if we use ranking method as given in 
Definition 3.2 and Definition 3.6 to compare two INNs, 
the results might be irrational. We can use the order 
relation <n to compare two INNs if the comparing 
results is obviously irrational by using the order rela- 
tion <p and vice versa. Therefore, the order relation 
<H and <p can be viewed as mutually complementary. 
From Example 3.3, it is easy to obtain that the order 
relation <g is order-preserving for operations ranking 
method. Meantime, it can avoid the irrational results as 
mentioned in Example 3.1 and Example 3.2, therefore, 
comparing with <h and <p, the order relation <g is 
more superior. 



neutrosophic Choquet integral 
operator and its properties 

4.1. Fuzzy measure and Choquet integral 

Definition 4.1. [9, 12, 13, 20, 23] A fuzzy measure 
on X is a set function p : P(X) -> [0, 1] satisfies the 
following conditions. 

(1) /x(0) = 0, p{X) — 1 (boundary condition) 

(2) If A, B e P{X) and A C B then jl( A ) < p(B) 
(monotonicity) 

One can see that a fuzzy measure is a normal mono- 
tone set function which vanishes at the empty set. 
Furthermore, A fuzzy measure on X is said to be 

• additive if p(A U B) — ji{ A ) + p(B) for all dis- 
joint subsets A, B C. X. 

• cardinality-based if, for any A C X, pi A ) depends 
only on the cardinality of A. 

In the framework of the MCDM, X can be interpreted 
as a finite decision criteria set. p(A) can be viewed as 
the grade of subjective importance of subset ACI 
The monotonicity of the fuzzy measure means that the 
importance of a subset of criteria cannot decrease when 
new criteria are added to it [12], 
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Definition 4.2. Let f be a real-valued function on X, the 
Choquet integral of f with respect to a fuzzy measure 
p on X is defined as 

n 

(C) I fdp = y^Q(X(/)) - p(Xq + i )))/(*(,■)) 

i= 1 

where (•) indicates a permutation on X such that 
f(X( 1)) < / 0 ( 2 )) • < f{X(n) ) and X ( .) = {.*(,■), ■■■ , 

X(n)}> 2f(ii+l) — 0- 

When using a fuzzy measure to model the importance 
of the subsets of criteria, the Choquet integral can be 
viewed as a aggregation function [4, 24, 25]. It has been 
proposed by many authors as an adequate substitute to 
weighted arithmetic mean (WAM) or ordered weighted 
averaging (OWA) operator to aggregate interacting cri- 
teria [10, 14, 19]. The Choquet integral identifies with 
a WAM (resp. OWA) as soon as the fuzzy measure is 
additive (resp. cardinality-based). 

4.2. Interval neutrosophic numbers Choquet 
integral operator 

Definition 4.3. Let f : X — »• L be an interval neutro 
sophic number function on X. and p be a fuzzy measure 
on X. The interval neutrosophic number Choquet inte- 
gral {IN NCI) of f with respective to p is defined as 

n ^ 



on X. Then their aggregated value by using the IN NCI 
operator is also an interval neutrosophic number, and 



(C) / fdu = 



1 - 



n;u ( 



l - tI 

/(*«)) 



\ MX(,))-/x(X( J+ i)) 



1 - 



n;=. ( 



i - r. 



n; 



( jL \/i(X (l ))-^(X( J+ i)) 

1 Mif 



/(*«)) 

n 






( jU \^(X(i))-/4X(/ + i)) 



.(6) 



li: 






Ml 

where (•) indicates a permutation on X according to 
a given total order < such that /(x(i)) < f{X( 2 )) • • • < 
~f{x( n) ) and X { .) = {x (i -), • • • , X( n) }, X ( „ + d = 0. 

Proof. The first result can be directly obtained from 
Definition 4.3 and Theorem 2. 1 . Eq. (6) is easily proved 
by using mathematical induction on n. 



su, 

& 



Definition 4.4. Let f and g be two interval neutro- 
sophic number functions on X. f and g are said to be 
comonotonic about a given order relation < if 



(xi) < f(xj) iff g(xj) < g(xj) V/, j e {1, 2, 



(C) J f dp = ^2 [^(^(o) - M(^(i+t))] f(.xq)) (5) 

i= 1 



where (•) indicates a permutation on X according to 
a given total order < such that f(xp)) < /(x( 2 )) • • • < 
f(x ( „)) and X { .) — [xp), ■■■ , xpp), X(„ + d = 0. 

Example 4.1. Let X = {x\,X 2 }, an interval 
neutrosophic number function f : X —>■ L is given 
by f(x i)= ([0.59,0.68], [0.05,0.15], [0.10,0.18]), 
f(x 2 ) = ([0.65, 0.70], [0.25, 0.35], [0.15, 0.30]}. 
Since H{f{x i)) = 0.5682, H(f(x 2 )) = 0.5686, then 
f(xi) <H fix 2 ). 

Assume that fuzzy measure p 
given by p{0) — 0, /u.([jri }) = 0.5, 
p({x 1 , X 2 }) = 1. By Eq. (5), we have 

(C) f fdp 



X -» [0, 1] is 
li({x 2 }) = 0.3, 



The following propositions show some properties of 
the INNCI operator. 

Proposition 4.1. (Idempotency). Let f, g be interval 
neutrosophic number functions on X, and a be an INN. 
If fix i) — aforxj e X, then (C) / fdp — a. 

Proposition 4.2. Let f, g be interval neutrosophic num- 
ber functions on X and p, v be fuzzy measures on X. 
then 

(C) J fdp © (C) J fdv = (C) J fdip + v). 

Proposition 4.3. Let f be interval neutrosophic number 
function on X, a > 0, and p be fuzzy measures on X. 
then 

C) f fd(ap) = a(C) f fdp. 



= [/4(X) - p{X( 2 ))] 7(X(1)) + [p([x 2 }) - /4(0)] 7(X(2)) 

= ([0.6090, 0.6861], [0.0810, 0.1934], [0.1129, 0.2098]) 

Theorem 4.1. Let f : X — > L be an inten’al neutro- 
sophic number function on X, and p be a fuzzy measure 



Remark 4.1. According to Definition 4.3 and statement 
(4) of Proposition 2.1, it is easy to obtain Proposition 
4.1, Proposition 4.2 and Proposition 4.3. The set func- 
tions p + v and ap, a > 0, (in the Proposition 4.2 and 
4.3J, are not normal fuzzy measures ( see Definition 4.1) 
since their ranges are the inten’al [0, 2] and [0, a]. 
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However, we can also calculate the two expressions, 
(C) f fdp © (C) J fdv and (C) J fd(ap), by Eqs. (5) 
and (6). 

Proposition 4.4. (monotonicity). Let f, g be interval 
neutrosophic number functions on X. Fora St implied 
operation-invariant total order <, if f(xi)<g(xi) for all 
Xj G X, then 



for 1 < j < it), and tc ( 2> (J) — n(j) for 1 < j < io and 
70 + 1 < j < n. 

If ^K**7r( 2 )(7o+l))— ^('Lr( 2 t(7o+2))— " " " — ))’ 

then the function f (]> and g are comonotonic about 
the permutation tx (2) — ■■■ , 7 r (2) (n)). We can 

obtain (C) J f (2) dp<(C) f gdp . 

By the definition of 7 (2) , 



■ , Tt{n)) be the per- 
such that /(Xjr(i)) 



(C) / fdp<{C) / gdp. 



Proof. Let (jt(1), tt(2), ■ ■ 
mutations on (1,2, •••,«) 

7 : 2 )) 7 : ' ' ' fHf(Xn(n)\ 

If f and g are comonotonic about the permutations 
(7r(l), 7r(2), • • • , Ttiti)), then the conclusion is com- 
pleted. 

In the following, we consider the case that If f 
and g are not comonotonic about the permutation 
(7r(l), jt(2), • • • , 7t(n)). Let 



g(Xiz(i 0 )) = min {gixj.a))} 

1 <i<n 



1 < io < n 



The interval neutrosophic numbers function 7 (1) is 
defined as: 



/ ( ^-LrO')) 



( f(X7T(i 0 )), j < io, 

\ f(xn(j)). Otherwise. 



Let 7r (1) = (7r fl, (l), • • • , 7t (l) (n)) is a permutation 
such that |(x jr (i)(t))<|(x jr (i)( 2 ))< • • • <|(^d) (! - A for 
1 < j < io, and n (v Hj) = n(j) for ;'o + 1 < j < n. 

We will consider two cases. 

Case 1. If Mx n <O(i 0 +if<g(x n <O( i0+ ^^^% 
(x n a)^), then the function f ( 1 1 and g are comonotonic 
about the permutation = (7r (1) (l), • • • , TX (l \n)). 
We can obtain (C) f f (l) dp<(C) J gdp. 

By the definition of f (l \ 

(C) j fdp<(C) I f (l) dp<( C) I gdp. 

Case 2. If there exists at least one point Xj 0 such 
that g(x^(i) (; - 0+1) )>|(x jr (i) O - 0+1) ), we define the interval 
neutrosophic numbers function f (2> as 



(C) j fdp<(C) j f m dp< J f (2) dp<{C) J ~gdp. 
If there exists at least one point Xk 0 such that 

gCLr(2>(yto+l)) >§(* 7 r< 2 >(£o+l))’ re P eat the aboVe Ste P S 

until ko — n . Then we will define a series interval neu- 
trosophic numbers functions f (k> (k < n ) and permuta- 
tions 7t (k \ then the function and g are comonotonic 
about the permutation 7t (k> = (n (k H\), ■ ■ ■ , 

We can obtain (C) / f (k ^dp<(C) J gdp. 

By the definition of f <k K 

(C) J fdns<C) j f"’dii<- ■■<((-) 

J f (k) dp<{C) J gdp. 




Proposition 4.5. Let f, g be interval neutrosophic 
number functions on X. For a <h implied operation- 
invariant total order <, if f and g are comonotonic 
about order relation <, then 

(C) J (/ © ~ g )dn = (C) j fdp © (C) J gdp. 

Proof. According to Definition 4.4 and Definition 
3.4, it is easy to get the conclusion. 

Proposition 4.6. Let f be internal neutrosophic number 
function on X. For a St implied operation-invariant 
total order < and X > 0, 

(C) J Xfdp = X(C) J J dp. 

Proof. According to Definition 3.4, it is easy to get 
the conclusion. 

From Proposition 4.5 and 4.6, we can obtain the 
following corollary. 



~X 2 ) ( , _ J /(**( 7o))’ *0 + 1 < j < jo, 

n(J) \ J m (XnU)), Otherwise. 

Let 7 r (2) = ( 7 r ( 2 , (l), • • • , tt ( 2 ) 00) is a permutation 
such that |(x 7r ( 2 ) (( - 0+1) )<|(x 7r ( 2 )(( - 0+2) )< • • • <~g{x n a\ k) ) 



Corollary 4.1. Let f, g be interval neutrosophic num- 
ber functions on X, a be an INNs. For a <l implied 
operation-invariant total order < and X > 0, 



(C) 



X(f 0 a)dp — X 




fdp 0 d 
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It is not hard to see that the aggregation properties of 
the INNCI greatly depend on the given order relation. 

Based on the order relations <h and </>, the INNCI 
only has the Properties 4.1— 4.2. However, based on a 
<L implied operation-invariant total order, the INNCI 
has all the properties presented above, 

Zhang et al. [34] have developed interval neutro- 
sophic number weighted averaging (INNWA) operator 
and interval neutrosophic number weighted geometric 
(INNWG) operator, which are respectively defined as 
follows. 

Definition 4.5. Letaj (j — 1,2 ,n) be a collection 

ofINNs, and let INNWA : INN n INN, INNWA is 
called an interval neutrosophic number weighted aver- 
aging operator if 

n ® 

INNWA M (ai, ci 2 , ■ ■ ■ , a n ) = &>,a, 

i= 1 

where a> — (co i, co 2 , ••• ,(o n ) is the weight vector of 
di (i = 1, 2, • • • , n), with coi > 0 (i — 1, 2, • • • , n) and 

S“' = L Q 

Definition 4.6. Letaj (j = 1, 2, • • • , n) be a collection 
ofINNs, and let INNWG : INN" -+ INN, INNWG is 
called an interval neutrosophic number weighted geo- 
metric operator if 

n ® 

INNWG m {a \ , d 2 , ■ ■ ■ , ci n ) = ci\ f 

i= 1 

where a> = (co 1 , ( 02 , ■ ■ ■ , co n ) is the weight vector of 

a, (i = 1, 2, • • • , n), with co, > 0 (i = 1,^, •••,«) and 
n 

J2 COi = 1- _ 

i=l 

The following theorems show that INNCI operator 
is a generalization of INNWA and INNWG. 



<o n -k — p({i + 1, i + 2, • • • , i + k}) — for 

any {/ + 1, i + 2, • • • , i + k} C X, 

CO 1 = 1 — (W2, 

then the proof is completed. 

Remark 4.2. The relationship between fuzzy measure 
and weight vector can be given as follows: 

(1) The additive fuzzy measure /x associated to an 
INNWA is given by pt{A) — "YhxinA c °if or A C. X. 

(2) The cardinality-based fuzzy measure p associ- 
ated to an INNWG is given by p(A) — Y^i=n-\A\+\ Wi 
for any non-empty subset A C. X. 

5. Multieriteria decision-making method based 
on interval neutrosophic numbers Choquet 
integral operator 

This section presents a new method for MCDM, in 
which the partial evaluations of the alternatives are 
given by INNs and the interaction among the criteria 
are allowed. 

For a MCDM problem, let Y = { y 1 , >' 2 . • ■ ■ , y m } be a 
set of alternatives, X — [xi , X 2 , ■ ■ ■ , x„] be a set of cri- 
teria. To get the best alternative, the MCDM procedure 
based on INNCI operator is proposed as follows. 

Step 1. Construct the INNs decision matrix. Assume 
that the partial evaluation of the alternative yj (i — 
1 , 2, • • • , m) on the criteria xj (j — 1 , 2, • • • , n) is mea- 
sured by INN djj = ( Tp , Ip, Fp), where Tp indicates 
the degree to which the alternative yj satisfies the crite- 
rion Xj, fj indicates the indeterminacy degree to which 
the alternative yy satisfies or does not satisfy the crite- 
rion Xj, Fij indicates the degree to which the alternative 
v; does not satisfy the criterion xj, and Tp C [0, 1], 
Ip C [0, 1 ],Fp c [0, 1 ]. Then we can obtain a decision 
making matrix as follow: 



Theorem 4.2. Let f be an interval neutrosophic number 
function on X, p be a fuzzy measure on X. 

( 1 ) If p is additive, then there exists co e [0, 1 ]" , such 
that (C) f fdp = INNWA. 

(2) If p is cardinality-based, then there exists co e 
[0, 1]", such that (C) J fdp = INNWG. 

Proof. (1) If >u is additive, let co, — p({i\), then the 
conclusion is obvious. 

(2) If p is cardinality-based, co is defined as follows: 
co n — p({i}) for any {/'} C X, 

(o n - 1 = p({i, j }) - (o n for any {/, j} c X, 



/ d\\, du, ■ ■ ■ , d\ n \ 

d-21, (I22, ■■■ ,d2n 



\ d m \ , d m 2 , • • • , d mn ) 

Step 2. Reorder the partial evaluation dp of the 
alternative y j such that d^j) < d^j+i) for a given order 
relation. For example, by order relation <// of Defini- 
tion 3.2, we calculate H(dp) of the partial evaluation d,j 
of the alternative yj (j — 1 , 2, • • • , m) on the criteria x, 
(i = 1, 2, • • • , n ), and then rank the partial evaluation 
dp such that dpj) < H d i( j+\)- 
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Step 3. Identify the fuzzy measure on the criterion 
set X = {x i , X2 , ■ ■ ■ , x„ } . There are several methods for 
the determination of the fuzzy measure, such as linear 
methods [18], quadratic methods [6, 7], heuristic-based 
methods [8], genetic algorithms [28] and so on. 

Step 4. Choose the INNCI operator to aggregate the 
partial evaluations of each alternative and get the overall 
evaluations. 

Step 5. Rank those overall evaluations according to 
the given total order relation on INNs, and select the 
best one(s). 



6. An illustrative example 

This section presents an application of the proposed 
method to select the third party logistics providers. Sup- 
pose that there are four providers (vi, yi, yi, >’4) whose 
core competencies are evaluated by means of the fol- 
lowing four criteria (.rj , x%, X3, *4): 

(1) the cost of service (xi); 

(2) the operational experience in the industry (xx)', 

(3) customer satisfaction (X3); 

(4) market reputation (*4). 

6.1. Procedures of decision making based on three 
order relation <h> <p and <5 

Step 1. The INNs decision matrix of the third party 
logistics providers is made up according to the four 
evaluating criteria. The evaluation of a provider yfi = 
1, 2, 3, 4) with respect to a criteria xj(j = 1, 2, 3, 4) 
is obtained from the experts. Suppose that the INNs 
decision matrix is constructed as shown in Table 1. 

Step 2. According to Table 1 , by Definition 3 . 2, Defi- 
nition 3.5 and Definition 3.6, the partial evaluation of 
the candidate yj is reordered such that d^j) <* dj(j+ 1) 
(i — 1, 2, 3, 4), where <* denotes the order relations 
<H, <p and <s, respectively. The results are shown in 
Table 2. 

Therefore 

^i(i) = d 12, di(2) — dn, d 1(3) = dn,d\(4) — dw, 

d 2 (\) = djy, d2(2) = <?21> ^2(3) = dl3 , <?4(4) = d 22 ', 

^ 3 ( 1 ) = dn, d 2 (2) — dn, ^ 3 ( 3 ) = d 2 3, dy(4) — dn', 

d 4(1) = dn, 4/4(2) = dn, 4/4(3) = dn, 4/4(4) = 4/41 . 

Step 3. Suppose that the fuzzy measures of criteria 
of X and subsets of X are shown in Table 3. 



Table 1 

The INNs decision matrix of the third party logistics providers 



*1 


yi 


([0.60, 0.76], [0.22. 0.33], [0.23, 0.30]> 


,V2 


([0.66, 0.71], [0.26, 0.31], [0.25, 0.33]} 


y3 


([0.66, 0.71], [0.26, 0.31], [0.25, 0.33]} 


V4 


([0.60, 0.76], [0.22. 0.33], [0.23, 0.30]} 


x 2 


yi 


([0.70, 0.72], [0.31. 0.37], [0.28, 0.31]} 


,V2 


([0.68, 0.75], [0.19, 0.29], [0.28, 0.35]} 


.V3 


([0.68, 0.75], [0.19, 0.29], [0.28, 0.35]} 


V4 


([0.70, 0.72], [0.31. 0.37], [0.28, 0.31]} 


x 3 


yi 


([0.71, 0.73], [0.21, 0.27], [0.35, 0.39]} 


V2 


([0.55, 0.73], [0.20. 0.25], [0.23, 0.30]} 


yi 


([0.55, 0.73], [0.20, 0.25], [0.23, 0.30]} 


V4 


<[0.71, 0.73], [0.21, 0.27], [0.35, 0.39]} 


a: 4 


y i 


<[0.71, 0.77], [0.23, 0.30], [0.29, 0.37]} 


yi 


([0.65, 0.70], [0.31, 0.32], [0.25, 0.27]} 


yi 


([0.65, 0.72], [0.25, 0.32], [0.22, 0.35]> 


y4 


([0.61, 0.74], [0.28, 0.30], [0.20, 0.32]} 



Table 2 

Ranking order of dq based on the order relation <*. 



ranking order 


<* 


Vl 


4l2 <* ^13 <* 4ll <* 3,4 


yi 


^24 <* ^21 <* ^23 <* d22 


yi 


dn <» 5 3 4 ^33 <* ^32 


V4 


^42 <* ^43 <* du <* d4i 



Table 3 

Every subset A of X and its corresponding fuzzy measure value 



A 


4(A) 


A 


M(A) 


A 


U(A) 


0 


0 


{3} 


0.10 


{1,3} 


0.20 


{1} 


0.10 


{4} 


0.10 


{1,4} 


0.50 


{2} 


0.10 


{1.2} 


0.30 


{2,3} 


0.4 


A 


4(A) 


A 


M(A) 






{2,4} 


0.40 


{1,3,4} 


0.60 






{3,4} 


0.20 


{2,3,4} 


0.30 






{1,2,3} 


0.85 


{1.2, 3, 4} 


1 







Table 4 

Results of aggregation based on INNCI operator 
Results of aggregation 

yi ([0.6657, 0.7428], [0.2523, 0.3354], [0.2656, 0.3187]) 

y 2 ([0.6308, 0.7189], [0.2391, 0.2901], [0.2466, 0.3130]) 

y 3 ([0.6335, 0.7193], [0.2338, 0.2878], [0.2498, 0.3207]) 

y 4 ([0.6582, 0.7362], [0.2461, 0.2965], [0.2624, 0.3430]) 



Step 4. By Eq. (6), utilizing the INNCI operator to 
aggregate the partial evaluations of each alternative y, 
(i — 1, 2, 3, 4). The results are shown in Table 4. 

Step 5. According to the overall evaluations of 
the third party logistic providers, by Definition 3.2, 
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Table 5 

Ranking order of the third party logistics providers. 



Table 6 

The unreasonable results based on the order relation <# and </>. 



Order relation 


Ranking order 




Corresponding unreasonable ranking order 


<H 


yi >- y 2 >- y3 >- j4 


yi >z y4 


yi <z y4 


<P 


yi >- yi >- y4 >- yi 


yi <z J 3 


yi >z yi 


<5 


yi >■ y4 >- ys >- J2 


yi <h ys 


yi >h yi 



Definition 3.5 and Definition 3.6, the ranking order 
based on the order relation <* of providers is shown 
in Table 5. 

6.2. Comparison analysis of results 

In Section 6. 1 , we have given the ranking order based 
on three order relations <h, <p and <$, but from the 
results shown in Table 5, the raking order calculated 
based on three order relations is different. Therefore, it 
is very difficult to decide which alternative is the best 
choice and whether an alternative is definitely better 
than another one. 

From Table 1 , we have that 



dll — d4 \ , d 12 = t/42, «13 = t/43. 



d 14 >S d-AA, 



d i4 > h d 44 (since H(di 4 ) — 0.7496 > 77(^44) 
= 0.7460), 



d \4 >p d .44 (since p{s(di 4 ) > 5(^44)) = 0.5417). 

Which mean that, on each criterion, the partial eval- 
uation of yi is larger or at least equal to that of V4 
(denoted as yj >- y4). Hence, the overall evaluation of 
y 1 should be larger or at least equal to that of 54. But 
unfortunately, aggregation function INNCI produce the 
contrary results based on the order relation <p : 

INNCI(yi ) < P INNCI(y 4 ). 

It is obvious that the above result is unreasonable. 
Table 6 shows the unreasonable results in this illus- 
trative example based on the order relation <h and 
<p. The main reason for such unreasonable results is 
that, as mentioned in Section 3.3, the order relation 
<H and <p are not order-preserving for operations, 
that is , these two order relations are not an operation- 
invariant total order on INNs. The increasingness 
property of an aggregation function strongly depends 
on the order-preserving property of the given order 
relation. Therefore, the interval neutrosophic numbers 



o 



Choquet integral operator, INNCI, is not increasing 
with respect to these two order relation. 

The key to avoiding such unreasonable results is to 
adopt a <l implied operation-invariant total order (see 
Definition 3.5) in the aggregation process. Based on the 
<L implied operation-invariant total order <5, we can 
generate the overall evaluations of the alternatives by 
using the INNCI, as shown in Table 5, and then we will 
select provider yi . 



7. Conclusions 

Interval neutrosophic set (INS) is a subclass of 
a neutrosophic set, which can be applied in the 
problems with uncertain, imprecise, incomplete, and 
inconsistent information existing in real applications. 

s an aggregation function, the Choquet integral 
with respect to fuzzy measures is able to flexibly 
describe the relative importance of decision criteria 
as well as their interactions. In this paper, we com- 
bined the Choquet integral and the INS theory to 
propose INNCI operator for MCDM problem with 
netrosophic information and investigated their aggrega- 
tion properties, such as idempotency and monotonicity. 
INNCI operator can represent INNWA and INNWG. 
Therefore, INNCI operator is superior to existing 
operators. 

Increasingness is a natural requirement for an aggre- 
gation function in MCDM. The increasingness property 
of INNCI aggregation function strongly depends on 
the order-preserving property of the given order rela- 
tion. In this paper, we proposed two approaches to 
compare two INNs. According to its geometrical struc- 
ture, the ranking index is developed according to 
some geometrical quantities which reflect the infe- 
riors and superiors of INNs. Based on the ranking 
index, an order relation, denoted by <#, is proposed. 
Examples shown that order relations <h and <p 
have some limitations, but they can be viewed as 
mutually complementary if the comparing results is 
obviously irrational. Because order relation <h and 
<P are not order-preserving for operations, we fur- 
thermore proposed a <p implied operation-invariant 
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total order to ensure the increasingness of INNCI oper- 
ator, and it is more superior than the order relations 
<H and <p. 

We only have proposed a kind of a </ implied 
operation-invariant total order (Definition 3.5), which 
does not adequately reflect the characteristics of INNs. 
Therefore, it is of great interest to find a <p implied 
operation-invariant total order which can better reflect 
the characteristics of INNs. Furthermore, how to find 
a good aggregation operation is also an important key 
issue in netrosophic MCDM problems. 
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